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We consider the bound states of the massive scalar field around a rotating black hole 



> 

o 

' immersed in the asymptotically uniform magnetic field. In the regime of slow black hole 

rotation, the Klein-Gordon equation allows separation of variables. We show that the 
growth rate of the instability can be amplified a few times by the magnetic field. The 
effect occurs because the magnetic field adds the "effective mass" term B\m\ to the scalar 
field potential for a Kerr black hole. In addition, and as a by-product, we discuss the 
behavior of the quasinormal modes for the magnetized rotating black holes. 
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fi; 1 Introduction 

■ Interaction of black holes and magnetic fields is an important factor for large astrophysical 
^ . black holes, and, especially, for supermassive galactic black holes because of enormous magnetic 
^ ! fields in the nucleus of galaxies and near supermassive black holes. Strong magnetic fields are 
O I induced also by accretion disks of rotating charged matter near black holes. Yet, interaction of 
strong magnetic fields and black holes is, apparently, not limited just by the context of large 
scale astrophysical systems. Recent development in the brane world theories opens possibility 
O ■ of observing quantum gravity at the Tev scale, so that mini black holes might be created at 
particle collisions in the Large Hadron Collider. Observing black holes in laboratories could 
make it possible to test interaction of black holes with external fields. After all, in the early 



^ ■ universe strong magnetic fields could effect primordial black holes. 



Sufficiently strong magnetic field near a black hole deforms the black hole space-time, so 
that one cannot consider the magnetic field as a test field on the black hole background, but a 
consistent solution of the Einstein-Maxwell equations is necessary instead. Fortunately, exact 
solutions for such a situation exist for both non-rotating (the Ernst solution) [1] and rotating 
cases (the Diaz solution) [2]. 

Yet, the essential question is whether such a strong magnetic field, that could deform the 
black hole space-time, can exist, and, whether the Ernst and Diaz solutions can be considered 
as at least some approximate approach to a realistic situation?. The Ernst solution is described 
by the metric 

,2 ,off 2M\,2 A 2M\"\ 2 2 ,n2^ r^sin^^,. 
ds' = A'Ul-—jdt'-i^l-—j dr'-r'de'j-^^d<j>^ (1) 

where the external magnetic field is determined by the parameter B, 

A = 1 + ^5^2 sin^ e, (2) 



and the "unit" magnetic field measured in Gs is 



Bm = 1/M = 2.4 X 1019^^. (3) 

The magnetic field in the Ernst solution is poloidal and homogeneous far from the black 
hole. This is not a big restriction for our consideration. Indeed, although the large scale 
magnetic field in the observable universe has both toroidal and poloidal components, it is 
the poloidal component that dominate in centers of galaxies and thereby in the region near a 
super-massive black hole. For mini black holes, which according to the brane-world scenarios, 
are expected to be observed in particles collisions, the homogeneous (poloidal) magnetic field 
is natural. Even though the super-strong accelerating magnetic field is screened in the region 
of particles collisions, we can assume, at least in principle, that if mini black holes are created 
in a laboratory, one could "immerse" mini-black holes in the magnetic field with required 
properties. 

In order to make estimations of possible infiuence of the magnetic field on the supermassive 
black holes we need the two parameters at hand: the magnetic field parameter B and the mass 
of the black hole M. Modern observations suggest that super-massive black holes in centers 
of galaxies can have mass M ^ 10^ — lO^Msun, while observations of the magnetic field are 
much more complicated and lead sometimes to controversial predictions for the value of the 
magnetic field in the centers of galaxies [3], At the same time, observations of the Active 
Galactic Nuclei (AGN) and micro- quasars indicate the existence of wide X-ray emission of 
lines of heavy ionized elements in their spectra [3j. This effect is very well described when 
supposing existence of the super-strong magnetic field B 10^° — 10^^ Gs |3]. In addition 
some other theories [S] imply existence of a strong magnetic field in the AGN. From these data 
for M and B at hand, by formula (2), one can easily see that for the heaviest galactic black 
holes, the magnetic field can be of order of fractions of a "unity" (1/M), thereby effecting the 
black hole metric itself. 

The Diaz and Ernst solutions, which we shall analyze here, describe the black holes im- 
mersed in a magnetic field decaying near the black hole to some asymptotic value, so that far 
from black hole the magnetic field is uniform and represents the Melvin Universe. Different 
effects around such black holes have been studied in [6] and some generalizations were obtained 
in [7]. Potentially observable effects, such as quasinormal modes and gravitational lensing were 
considered in [8], [9], [10] for non-rotating magnetized black holes. In this paper, we shall con- 
sider scattering of scalar waves around rotating magnetized black holes. An incident wave near 
the black hole will be partially absorbed (tunneled through the potential barrier) and partially 
radiated away as a response of a black hole to an external perturbation. When a black hole of 
mass M and radius r+ is rotating, and the real oscillation frequency of the perturbation (for a 
mode with an azimuthal number m) satisfies the inequality 

am , 
u < , (4) 

then the energy radiated away exceeds the energy of the incident wave. This is the well known 
effect of super-radiance predicted by Zeldovich and Misner [TI] and computed for Kerr black 
holes by Starobinsky [12]. Press and Teukolsky suggested that if the wave radiated away will 
be reflected by a mirror surrounding a black hole, one could make a kind of "black hole bomb", 
because initial small perturbations would grow without bound [13]. The role of the mirror can 
play the potential well, another local minimum, which appears because of the non-vanishing 
massive term 
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In the limit M/i << 1, where M is the black hole mass and jj is the inverse Compton wave- 
length of the particle, Detweiler showed that the massive scalar field exerts the superradiant 
instability with the maximal growth rate ^15j, ^6] (at a = M) given by the formula 

This instability has very small growth rate and probably cannot be significant for real black 
holes. Accurate calculations of instability made recently in [T7] allowed to observe the maximal 
growth 

7 = 1.510-^(GM/c3)-\ (6) 

at M/i = 0.42. Yet, due-to short lifetimes or large particle masses, this instability is tiny for 
any known massive boson particles, except possibly, the neutral pion 7r° when M ~ 10^^%. 

In this paper we shall consider the bound states of massive scalar field in the vicinity of 
the magnetized black holes and analyze the corresponding superradiant instability. As a by- 
product we shall consider the 5-corrections to the orbital geodesic motion (as that for which 
the instability is important) of particles around the Ernst black hole. The paper is organized 
as follows. In Sec. II the massive scalar field equation is decoupled for the regime of "small" 
(in comparison with a "unity" 1/M) magnetic field, and is reduced to the wave-like equation 
with an effective potential. Sec. Ill considers the bound states and their instability in the limit 
/iM -C 1. In particular, we obtain the generalization of the well-known Detweiler formula for 
the superradiant instability for a non-zero external magnetic field. In Sec. IV we briefiy discuss 
the quasinormal modes of rotating magnetized black holes. In the Conclusion we summarize 
the obtained results and outline the number of open questions. 



2 The metric and separation of variables in the wave 
equation. 

The Diaz metric is the exact solution of the Einstein-Maxwell equation and represents a charged 
rotating Kerr- Newman black holes immersed in an asymptotically uniform electromagnetic field 
(se formulas (2-11) in [18] and [2]). If we put the black hole charge, and external electric field 



to be zero, then the metric takes the form 

els'" = lAI^S (^jdt' - A-'dr^ - de'^ - ^^id(j) - Wdtf. (7) 

In order to separte variable in the Klein-Gordon equation, we have to be limited by the case 
of small magnetic fields B « M^^, and small rotation a « M, so that we neglect terms of 
order B'^a'^ or B^ and higher. Within this approximation 

A = l + ^5Vsin2e + 0(5V,5^), (8) 

W=^-^ + OiB'a^B% (9) 
A = r^ + - 2Mr, (10) 
S = + cos^ ^, (11) 
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A={r^ + a^f-/\a^siY?e. (12) 

Here B is the parameter of the magnetic field, which gives the asymptotic value of the 
magnetic field. The solution (7), as was shown by Aliev and Galtsov [19], forms the conical 
singularity on a polar axis. This generates additional (to the Maxwell term) singular term in 
the stress-energy tensor on the right hand side of the Einstein equations. In order to avoid this 
singular term, one should consider the more general metric given in [19] (see equation (2.9) in 
[19]). Yet in the slow rotation regime and zero background charge (5 = 0, the Aliev-Galtsov 
metric reduces to the uncharged Diaz metric, which we shall use here. 

The massive scalar field obeys the equation 

^ ^ (13) 



As the background metric has the Killing vectors {d^f), dt), it is implied that 

$ ~ e^™'^+*'^*5(^)i?(r). (14) 

The wave equation (13) consists of the five terms, corresponding to components of (13) with 
tt, rr, 66, t(j), and 00 derivatives. The dependence on A (and thereby on B) is preserved only 
in the last term 00, while in the other terms, A~^ coming from the contra- variant components 
of the metric is canceled by A^ of \f~-g- Neglecting small terms of orders 0{B'^a'^,B^,B'^fi'^) 
and higher, one can reduce the wave equation to separate equations for S{6) and R{r) in the 
limit of small rotation and "weak" magnetic fields 



A A ( A^V(^'(^' + «')'-4aMrmcu + aW-A((5W + /x2)r2 + aV + A))i? = 0, (15) 
dr \ dr J 

1 d /gij^^^MV (a2^2^^g2^_^2gi^-2^^^)^(^)^ ^^g) 



sin^rf^ V 

Thus we see that the radial equation has the form of wave equation for the Kerr black hole 
with the "effective mass" \/ ^ B'^w? , while the equation for angular part looks like that for 
the massless scalar field, because the "massive" term in the angular part is of orders a?B^. Yet, 
that is not a problem for our future analysis, because it is enough to approximate A as £(£ + 1). 
In the limit a = we have the effective potential of the non- rotating case pLOj. It is worthwhile 
to note that in a similar fashion, the effective mass term appear for the massless scalar field 
when considering dominant term of the self-action of the scalar field or when a black hole is of 
Kaluza-Klein type [20]. 

Let us stress here that the "small" magnetic field is a very good approximation to the real 
situation and means here in fact a super-strong magnetic field. Indeed, as the " unit" magnetic 
field measured in Gs is Bm = 1/M = 2.4 x 10-^^ (M5„„/M), and if we take M = 1 for, say, a 
black hole of mass l^'^Msuni then B = 1 means B = 2.4 x 10^° Gs, what is about the possible 
maximum one can expect in astrophysical context [3]. We made one implicit approximation, 
when discarding the terms of order 5^ and higher: in fact we discarded terms S^r'^ and higher, 
so that we can trust our equations (15, 16) only for Br <^ 1 (M = 1). This again is not a 
serious restriction as even for a super-strong field, S is a few orders less than unity and thus 
the equations (15, 16) are a good approximation for r ranging from the event horizon until 
1/B, what is few orders bigger than the radius of the event horizon. In this way one has 
"infinity" situated somewhere around r = 1/B and we solve thereby another problem: the 
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exact Ernst or Diaz metrics are not asymptotically flat, because the magnetic field exists even 
at asymptotic infinity. This is certainly not physical situation and the Ernst metric should 
have been matched somewhere far from the black hole to an asymptotically fiat solution. We 
do not need to do it here because the scattering processes happen near the peak of the effective 
potential barrier, and practically, incident wave comes to barrier not from asymptotic infinity 
but from the far region r ^ = 2M. Therefore the scattering properties, such as absorption 
probabilities or quasinormal modes are determined by the behavior of the effective potential 
near its maximum. In the astrophysical context, the independence of the quasinormal modes 
on the asymptotic behavior of the effective potential was proved in |21j. (This certainly does 
not concern asymptotically AdS space-time where the boundary conditions at spatial infinity 
is interpreted through the AdS/CFT correspondence). 



3 Bound states and superradiant instability 

The bound states appear when the effective potential has a local minimum. For M = 1 and 
different values of i, bound states are formed for all values of the effective mass term less than 
some threshold value \/ jJ.'^ + B'^w? < (v^yU^ + B'^m?)max- The values of (^/x^ + B'^w?)max for 
different values of ^ are plotted in Fig. 1. For large ^ this maximal threshold value is 

(V/i^ + ^W)^,, = v7(7+l)/2y3, (M = 1) (17) 



ngular number 1 



Figure 1: The maximal value of the effective mass term MyZ/iM-^B^m^ allowing the 
bound metastable states for different values of £. One can see that the maximal values of 
M-\//l2~+~B%i2 approaches its asymptotic value /2\/?)M in the regime of high ^ 

already at moderate £ 

From Fig. 1 one can see that this asymptotic regime is reached already at moderate I. The 
discussion of the possible values of particle and black hole parameters when bound states are 
formed was given in [22] for the Schwarzschild black hole. 

The instability of the bound states is not very important if a particle moves towards a 
black hole: the particle will simply fall into the black hole. Yet for particles moving around 
black holes on classically stable orbits, quantum metastability determines the time of orbiting 
around the black hole. Therefore in the Appendix we shall find B - corrections to the classical 
orbital motion of particles around Ernst black holes. 

Let us start from the boundary conditions for the bound states problem 

^ ~ 0, — > +00, 
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which are the pure in-going waves at the event horizon and zero at the spatial infinity. Under 
these boundary conditions the massive scalar field has bound states, which "live" in the local 
minimum of the effective potential barrier. 

Here we shall follow the Detweiler approach which uses the approximation /iM << 1. In 
our case, in addition to this approximation, there are the two other restrictions: small rotation 
a <^ 1, and "small" magnetic field B ^ 1, (if M = 1). We shall not repeat all the steps of the 
Detweiler work [15], but just say that they are the same if one change fi"^ by the /j,^ + B'^m'^. 
This leads to the following result for the frequency of the bound state solution of the above 
wave equation. 



uj = (a/ jj,"^ + B'^m 



+ 27, (18) 



where the imaginary part 7 describes the rate of growth (or of decreasing) of each mode, 

7 = + B^m^{M{^fi^ + B^m?)Y^^\amM-^ - fi^ + BV)r+)x 
2^^+2(2^ + 1 + n)! / a ' ^ 



X 



n!(£+ 1 + n)2^+4 \{2i)\{2e + ly. J 
n5=i(j'(l - a'/M^) + ((am/M) - 2{^/ fi^ + B^m^)r+f). (19) 
If one has m > and 



amM-^ - 2{^JFTB^)r+ > 0, 5 < Bm (20) 

then 7 > and we have a super-radiant instability. The analog of the 2p state of the hydrogen 
atom (£ = m = 1, = 0) has the fastest instability growth, if we are limited by the regime 
BM\m\ « 1. Due-to condition ^'^B'^ -C 1, in the obtained formulas (18,19) we are limited by 
those values of the magnetic field B and azimuthal number m, for which B\m\ at least a few 
times smaller than /x, so that the contribution due-to magnetic field plays a role of a relatively 
small (but not negligible) correction to the effective mass of the scalar field. In this case, the 
"asymptotic" region 1/i? is situated at a distance which is a few times larger than the distance 
from the black hole to the local minimum (see Fig. 2). 

The growth rate as a function of the magnetic field can be found in Fig. 3. There one can 
see that in the region of stability, the magnetic field increases the damping rate of a mode, 
while in the region of superradiant instability, magnetic field can increase the growth rate 
by a few times. Let us stress that this result was obtained for the slow rotation regime. It 
is well known that for highly rotating black holes the superradiant instability growth rate is 
increased by many orders [18], [15]. For highly rotating black holes the coupling between the 
projection of the angular momentum, given by m, and the magnetic field B, must be much 
stronger, leading, probably, to a stronger superradiant instability. Anyway, let us speculate 
what happens if at least a similar, few times increasing of the growth rate occurs for a highly 
rotating black holes. For a superradiant instability to be significant in real processes it is 
important that the e-folding time of the maximal instability growth be larger than each of the 
two quantities: the lifetime of the black hole and the inverse decay rate of unstable particles. 
The accurate estimations given in [17] for a highly rotating Kerr black holes show that the 
lifetime of neutral pion T('7ro) = 8 x 10~^^ s. is of the same order as the e-folding time for the 
maximal instability growth r = 1.5 x 10"^'' s. Yet this happen only in the very narrow region 
of masses of black holes M ~ 10^^%. Assuming few times increasing of the instability growth 
rate in the presence of the magnetic field, we could expect much stronger instability for neutral 
pion and much wider range of black hole masses for which such instability is essential. 



6 



20 



40 



60 



80 



0.09 - 



0.08 - 

Figure 2: The effective potential m = £ = 1, B = O.Of, fi = 0.3, M = 1. One can trust the 
effective potential in the the region (2, 1/B). The formal local minimum, where the bound 
states can hve, appears at r ~ 20M. 



Figure 3: The growth/damping rate 7 as a function of the dimensionless parameter BM for 
massive scalar field: i = m = l, n = 0, a = 0.1, fiM = 0.1 (left), /xM = 0.05, fxM = 0.01 
(right). One can see that there is no superradiant instability for cases = 0.1 and 0.05 as 
for those beyond the inequality (20). The superradiant instability growth rate, plotted in the 
right figure, is shown to be increased by a few times because of the presence of the magnetic 
field. 
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4 Quasinormal modes 



Let us discuss here the propagation of the both massive and massless scalar fields with another 
boundary conditions, called the quasinormalhoundaij conditions [23] ■ These are pure in-going 
waves at the event horizon and pure out-going at the spatial infinity. 

^^e"*'^''*, -oo, and ^ ~ e*'^'^*, +00, 

The absence of the waves incident on a potential barrier of the black hole gives rise to the 
formal absence of superradiance. Therefore one cannot see superradiant instability under pure 
quasinormal boundary conditions. After all, as it is clear from the previous section, the insta- 
bility is significant only for a narrow class of particles and black hole masses. The quasinormal 
modes of the considered here slowly rotating magnetized black holes reduce to the quasinormal 
modes of Kerr black holes with "re-defined" mass of the field /i 1^/1^ + B'^m?. It is essential 
that, similar to the scattering characteristics, the quasinormal modes also depend on the be- 
havior of the effective potential near its peak, while behavior of the potential far from the black 
hole is insignificant [21]. This was proved in [21] by showing that quasinormal modes of two 
different potentials, with equal behavior near the black hole and strikingly different behavior 
in the far region, coincide. Using this fact, the quasinormal modes of non-rotating magnetized 
black holes were found in [1^. In a similar fashion, one can re-obtain the quasinormal modes 
of the slowly rotating black holes in the magnetic field by using numerical data of [21]- [26] and 
by the re-definition yU A//i^ + B'^nrfi. Here we shall not repeat these stages but refer a reader 
to the work [24|. Here we simply enumerate the main features of the quasinormal spectrum: 
First, B\m\ acting like a "mass", leads to decreasing of the imaginary part of the quasinormal 
mode until. 

The greater i?|m|, the less the damping rate of the mode and for sufficiently large value of 
Y^/i^ + B'^rn?, the infinitely long lived modes, called quasi-resonances p5], may appear. Let 
us remind that we separated varibles under the condition /iM -C 1, so that the existence of 
the quasi-resonances (which happen when /xM >~ 1) are guaranteed only for massless scalar 
field in a magnetized black hole background. Similar massless scalar field quasi-resonances 
were found for non-rotating magnetized black holes in [TO] and originally, for the Schwarzschild 
black holes in [26] . Second, in the regime of infinitely high damping, the quasinormal modes 
do not depend on magnetic field and are the same as for the massless scalar field in the pure 
Kerr background. Finally, the real oscillation frequency Reuj decreases when increasing Bm 
until some minimal value and then increases until reaching its quasi-resonance value. 

Here we see that the perturbations under pure quasinormal boundary conditions are longer 
lived in the presence of the magnetic field. Unlike quasinormal modes, scattering boundary 
conditions lead to the increasing of the damping rate of each mode in the stable sector (see 
Fig. 2). We know that the superradiant instability has the tiny growth rate, except some 
special cases of the neutral pion decay, which we discussed in the previous section. All this 
does not mean any contradiction. First of all, because quasinormal modes are beyond the region 
of superradiant instability. In addition, we implied here that, putting aside those exceptional 
cases, such an instability (and more general, superradiance) does not infiuence the quasinormal 
modes. Although any realistic perturbation contains some remnants of the incident wave of 
the initial perturbative process, at sufficiently late time the effect of this remnant is negligible. 

5 Conclusion 

We have separated variables for the Klein-Gordon and Hamilton- Jacoby equations in the limit 
B <C . This corresponds to the "weak" magnetic field as for the deformation the black hole 
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space-time, but a super-strong field in the astrophysical context. It is shown that the growth 
rate of the superradiant instability may be amplified by a few times due-to the magnetic field. 
This enlarges the range of black hole masses for which the instability may be significant in the 
neutral pion decay. The essential problem which was beyond our investigation is the calculation 
of the instability growth for the highly rotating black holes. Another important question is a 
detailed description of the process of filling of the quasi-bound states 



Appendix: Classical motion of particles near the magne- 
tized black hole 

The equatorial motion of particles in the background of the Ernst space-time was considered 
in [27]. The Ernst matric can be obtained if one takes a = limit in the Diaz metric (7). The 
Hamilton- Jacoby equation takes the form: 

dS dS 2 

Implying that S = —Et + L(j) + S'i(r, 6), one has 

. A (dS,V 1 (dS,V A' ^2 ,2,2 



r^[dr ) r^[d9 ) rHmH^ " ^ ' ^^^^ 

Here we shall not be limited by equatorial plane but will imply that B <^ and ^^B"^ <C 1, 
so that terms of order /i^i?^ <^ 1 and higher will be discarded. Thus we replace on the left 
hand side of equation (22) by 1 + i?^r^ sin^ 9, and A^ by unity on the right hand side. 

Then the effective potential takes the same form as a that for the Schwarzschild black hole 
with shifter angulat momentum — L^(l + i?^), 

r/2 ^ M f 1 + £!(1±£!) V (23) 

e// ^2 ^2^2 J 

Our main interest here is to determine parameters region when bound states, that is, 
classically stable orbits are possible. One can easily see that minimal distance r = 6M at 
which the stable circular orbit exists is unchanged in the lowest dominant order which we 
considered. Yet the minimal value of the angular momentum now depends on B 



""'"^ = 3.464. 

M 

In the equatorial plane the Hamilton- Jacoby equation (22) allows separation of variable without 
any approximation p7] and one can obtain the exact B - correction to the position of the stable 
circular orbit. 
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